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ABSTRACT 

We present the results of a 5.5-year CCD photometric campaign that monitored 
261 bright, southern, semi-regular variables with relatively precise Hipparcos paral- 
laxes. The data are supplemented with independent photoelectric observations of 34 
of the brightest stars, including 11 that were not part of the CCD survey, and a previ- 
ously unpublished long time-series of VZ Cam. Pulsation periods and amplitudes are 
established for 247 of these stars, the majority of which have not been determined be- 
fore. All M giants with sufficient observations for period determination are found to be 
variable, with 87% of the sample (at S/N ^ 7.5) exhibiting multi-periodic behaviour. 
The period ratios of local SRVs are in excellent agreement with those in the Large 
Magellanic Cloud. Apparent if -band magnitudes are extracted from multiple NIR 
catalogues and analysed to determine the most reliable values. We review the effects 
of interstellar and circumstellar extinction and calculate absolute if-band magnitudes 
using revised Hipparcos parallaxes. 

Key words: stars: AGB and post-AGB - stars: fundamental parameters - stars: 
late-type - stars: mass-loss - stars: variables: other - solar neighbourhood 



1 INTRODUCTION 

Photometric surveys of M giants in the Large Magellanic 
Cloud (LMC), Small Magellanic Cloud (SMC), and Galactic 
Bulge have revealed a rich variety of pulsational behaviour. 
Pulsating M giants are potentially important as distance in- 
dicators since they obey period-luminosity (P-L) relations. 
The P-L relation for Miras has been known for decades 



cise Hipparcos parallaxes, Bedding & Zijlstra ( 1998 1 found 



I Eggen|1975l [Robertson fc Feast |198T| [Glass fc Lloyd-Evans 



1981). More recent studies have shown that semi-regular 



variables (SRVs) also follow well-defined P-L relations. They 
may pulsate on the Mira sequence or in higher overtones, 
where they exhibit shorter periods and smaller amplitudes 
than the Miras (see, for example, Wood fc Sebo|1996 [Bed- 
ding & Zijlstra 1998'; Wood et al. "1999': 'Kis s fc Bedding 
2003, 20 04: Ita ct al. 2004. Groenewcgcn ,2004| [Wray Eyer 
fc Paczyn ski 2004; [Soszyiiski et al.|2007| ). 

Meanwhile, studies of M giants in the local solar neigh- 
bourhood have been hampered by imprecise distances and 
a lack of detailed, long-term observations to enable the de- 
termination of characteristic periods. Using published peri- 
ods for 24 galactic Miras and SRVs having relatively pre- 
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two parallel sequences in the {K, \ogP) plane, similar to 
those in the LMC. Further work by Knapp et al. ( 2003 1 and 



Ye§ilyaprak fc Asian (20041 increased the sample size and 



found evidence for at least two sequences, but the scatter 
due to imprecise parallaxes made it difficult to assign stars 
to distinct sequences. In the most extensive study to date, 
[Glass fc van Leeuwen[ ( [2007| ) used revised Hipparcos paral- 
laxes and literature periods for 64 red giants, confirming that 
local stars obey similar P-L sequences to those found in the 
Magellanic Clouds and Bulge. However, they concluded that 
details remain sketchy due to insufficient data, particularly 
for short-period stars. 

Due to their complex multi-periodic behaviour, SRVs 
require monitoring over several years to determine their peri- 
ods, particularly since many exhibit long secondary periods 
superposed with short-period, small-amplitude variations. 
Unfortunately, targets with good Hipparcos parallaxes are 
too bright for surveys such as AS AS ( PojmaiiskI|[1997 1 and 
NSVS ( [Wozniak et al.|2004[ ). 

We have constructed a wide-field, automated survey 
telescope optimised for precise photometry of bright local 
SRVs, with the dual goals of determining periods for a suf- 
ficiently large sample to enable unambiguous identification 
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of the local sequences, and to determine the effect of metal- 
licity on the P-L zero-point. In this paper, we describe our 
sample selection (Sect. [2]|, the selection of K-hand magni- 
tudes (Sect.jsJ, CCD instrumentation (Sect.[4]|, photometric 
pipeline and supplementary information (Sect.[5|, and PEP 
photometry (Sect.[6|. The main result of the survey is pre- 
sented in Sect. [7j where we provide periods and if-band 
magnitudes, and discuss variability, multi-periodicity, and 
period ratios. We summarise our results in Section |8] A de- 
tailed analysis of the results will be described in a follow-up 
paper, hereafter Paper II. 



2 SAMPLE SELECTION 

Our sample of nearby variable M giants was extracted from 
the Hipparcos catalogue (Perryman et al. 19971, using the 
following criteria: 



1997 



(i) a,r/7r ^ 0.2. Catalogue entries with negative paral- 
laxes were ignored. 

(ii) Coarse variability flag (H6) 2, indicating that the 
Hipparcos mission had observed amplitudes ^ 0.06 mag. 
This was used to select true variables. 

(iii) Spectral type (H76) = 'M*'. 

This yielded a complete sample of variable M giants to 20% 
parallax precision. A subsequent check prior to publication 
revealed that two stars (HIP 56518 & 85302) had been er- 
roneously omitted from the sample. 

The candidates were further filtered to select those 
likely to yield high-quality photometry: 

(i) Stars with declination > 35° were removed to ensure 
that all candidates would transit at altitude ^ 20° at the 
observing site. 

(ii) Stars fainter than mv = 8.0 at maximum light were 
excluded, retaining only those bright enough to yield pho- 
tometry at the 1% level. 

(iii) 33 stars listed in the General Catalogue of Variable 
Stars (GCVS) as having established periods > 100 d were 
removed in order to concentrate on those undergoing short- 
period, low-amplitude pulsation. 

This selection yielded 183 candidates. A further 67 stars 
with o",r/7r < 0.26 were added to increase the total number 
to 250, to provide sufficient targets for a 5-hour nightly ob- 
serving run. Stars were evenly distributed in right ascension 
so that a similar number of fields would be observed each 
night, regardless of the season. 

In hindsight, the declination-based selection criterion 
was too loose, resulting in the selection of a few stars with 
northern declinations that were only visible for short inter- 
vals each year and so have insufficient observations to de- 
termine their periods. In addition, one star near the south 
celestial pole (NSV 18694) was unobservable due to a phys- 
ical obstruction near the observatory. 



Since our original sample selection, van Leeuwen ( 2007 1 



has published revised Hipparcos parallaxes based on a re- 
reduction of the original data using an improved understand- 
ing of the rotation rates of the spacecraft (van Leeuwen 



& Fantino 2005). The version published on VizieR on 15 
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Figure 1. Distributions of relative uncertainty in parallax for the 
250 M giant sample, using the original (1997) and revised (2007) 
Hipparcos catalogues, binned into 5% intervals. 



sis in this paper. The original and revised distributions of 
the relative parallax uncertainties of our sample are shown 
in Figure [l] The revised distribution clearly shows an im- 
provement for the majority of stars, with the median error 
reduced from 17.5% to 10.6%, consistent with the finding of 



Glass & van Leeuwen (20071 that parallaxes for M giants 



were generally improved by a factor of 2. We have adopted 
the revised parallaxes for our analysis. 

Our CCD data are supplemented by the following 
sources: 

(i) Twenty-three stars were observed concurrently by one 
of us (T.T.M.) using photoelectric photometry (PEP), as 
described in Section [S] An additional 11 stars were observed 
exclusively by PEP, including 4 stars of spectral type C, one 
S, and one K-giant (A Vel). 

(ii) VZ Cam is a northern star with a relatively precise 
parallax measurement (4%). We include a previously unpub- 
lished light curve obtained by B.Sz. at Konkoly Observatory. 

Note that OO Aps was initially observed by PEP only, but 
was later added to the CCD target list to check for the 
possible existence of short-period pulsations. 



3 /^-BAND MAGNITUDES 

In preparation for a Period-Luminosity (P-L) analysis (Pa- 
per II), we obtained apparent K magnitudes as follows. 



3.1 2MASS 



September 2008, which corrected an error that affected ear- 
lier versions of the catalogue, has been used for the analy- 



2MASS ( Cutri et al.|2003 1 has been the catalogue of choice 
for previous studies of pulsating red giants in the LMC, 
SMC, and Bulge. However, most of our M giants were satu- 
rated and have large uncertainties of 0.2-0.3 mag. Moreover, 
2MASS magnitudes were derived from single-epoch mea- 
surements, leading to some variability-induced scatter, given 
that measured SRVs have -ft'-band peak-to-peak amplitudes 
between 0.1 and 0.25 mag ( [Smith et "ar| |2004[). Thus, we 



only used 2MASS magnitudes when better values could not 
be found. 
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3.2 The Gezari et al. Catalogue 

173 M giants in our sample appear in the Catalogue of In- 



frared Observations ( Gezari et al. 1999 1. All magnitudes and 



colour magnitudes (fiags M and C, respectively) observed at 
a wavelength of 2.2 ± O.OS^m were combined to calculate a 
mean rriK for each star, weighted equally because the cata- 
logue did not provide measurement uncertainties. 

Of the 80 stars with more than one catalogued mag- 
nitude, the mean scatter (excluding seven outliers) was 
0.04 ± 0.02 mag. Seven stars had an abnormally large scat- 
ter. In two, this was caused by an apparent transcription 
error where the sign of the magnitude had been reversed 
[jj. Gem and DM Eri). In other cases, a single discordant 
value, though not obviously wrong, was sufficiently different 
to yield a large scatter. As a precaution, mean magnitudes 
with a > 0.1 mag were not used. 



3.3 DIRBE 

During its 10-month cryogenic mission, DIRBE collected 
100-1900 measurements for each of nearly 12000 objects 



( Smith et al.|2004 l. This permitted an accurate determina- 
tion of mean K magnitudes for short-period SRVs {omK < 
0.06 mag). Furthermore, DIRBE was less prone to satura- 
tion than 2MASS. 

However, DIRBE had poor angular resolution due to its 
large beam width of 0.7°, which led to confusion between 
nearby NIR sources, particularly near the galactic plane 
(l&l < 5°). A search of the catalogue identified matches for 
222 project stars, but only 166 were marked as unconfused 
sources. In 5 cases, quoted errors exceeded that of 2MASS 
(due to insufficient observations) and were thus excluded. 

we adopt a value of 630 



Following Whitelock et al 
Jy for rriK 



|1988 


1. 


gure 





0, which is in good agreement with |Bessell fc| 



the mean rriK values derived from the Gezari catalogue, ex- 
cluding stars marked as confused in DIRBE and those stars 
extracted from Gezari with cr > 0.1 mag. The agreement is 
excellent, with a mean difference of —0.004 ±0.081 mag for 
the 125 stars on the plot, confirming that the calibration is 
appropriate. The small scatter is an indication of the consis- 
tency between the homogeneous DIRBE measurements and 
the mean values derived from multiple measurements taken 
from a set of inhomogeneous catalogues. Hence, it appears 
safe to use Gezari magnitudes when DIRBE values are un- 
available or confused with other sources. 

Figure |3] compares DIRBE magnitudes with 2MASS. 
The larger scatter than in Figure [2] is due to the saturation- 
induced errors in 2MASS magnitudes. The scatter for stars 
marked with confusion fiag 3 (001) is no worse than that 
for unconfused sources (000), so we retained both. However, 
stars with fiags 1 or 2 typically show an excess in the DIRBE 
measurement due to contamination with a nearby source, 
and have been excluded. 

In summary, DIRBE K-h&nA magnitudes were adopted 
except in cases where confusion fiag 1 or 2 was set, or 
the DIRBE uncertainty exceeded that of 2MASS. In those 
cases, Gezari mean magnitudes were used when available (41 
stars), otherwise the 2MASS magnitude was used (25 stars). 
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Figure 2. Apparent X-band magnitudes of unconfused DIRBE 
sources compared to mean Gezari magnitudes, excluding stars 
from the latter with rms scatter > 0.1 mag, showing an excellent 
correlation. 
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Figure 3. DIRBE vs. 2MASS X-band magnitudes using differ- 
ent symbols to denote various combinations of DIRBE confusion 
flags. The scatter is larger due to saturation in 2MASS. 



3.4 Interstellar Extinction 

Since all stars in the sample are nearby, with the most dis- 
tant being PS Hya (552 pc), we expect a relatively small 
amount of interstellar extinction, particularly in the K- 
band. Using & {B — V) vs. (V — 7) colour-colour diagram, 
we found that the majority of stars are unreddened. The 
maximum colour excess of E{B — 1/) ~ 0.2 occurred for 
7 stars with V — I > 3.5, implying a visual extinction of 
Av ~ 0.6 mag. This was confirmed using two independent 
interstellar extinction models ( Hakkila et al.|1997 and Drim- 



|mel, Cabrera-Lavers fc L6pez-Corredoira||2003[ ), which in- 
dicated a maximum of Av ~ 0.5 mag for all stars, with 
Av < 0.2 mag for the majority. Calculating if -band extinc- 
tion from Ak/Av = 0.112 ( |Rieke fc Lebofsky|1985,; Schleg^ 



et al. 1998), we find an upper limit of ~ 0.05 mag with 
Ak < 0.02 mag for 90% of the sample. Thus the correc- 
tion is small, particularly since Mk errors are dominated by 
much larger parallax uncertainties. Nevertheless, for the sake 
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also been used to discover a comet and two novae ( Tabur et 



al.|2002[[Tabur fc MonM^|2003a| [Tabur et al.|2003b i' 



Figure 4. Survey equipment: Nikon 80— 200-mm telephoto lens 
and SBIG ST8XE CCD detector attached to the telescope with 
a wooden cradle. 



of completeness and to avoid systematic errors, corrections 
were applied using the model of |Drimmel, Cabrera-Lavers| 
& Lopez-Corredoira (20031, due to its superior resolution. 



3.5 Circumstellar Extinction 

Small-amplitude SRVs have smaller circumstellar envelopes 
than their high-amplitude, more evolved counterparts, the 
Miras. To investigate the possibility of reddening due to 
a dusty envelope, we tested for IR excesses at 12/im by 
plotting V — I YS V — [12] for the 231 stars in the sam- 
ple having 12/im flux densities listed in the DIRBE cat- 
alogue. Only a few stars showed a small excess. A cata- 
logue of low-resolution spectra obtained by the Infrared As- 
tronomical Satellite (IRAS) confirmed that 8 stars had a 
small amount of circumstellar dust, while 15 showed silicate 



emission features. Using the relation of Knapp et al. ( 2003 



Equ 3) to estimate the -ff-band extinction caused by silicate 
dust, we found negligible amounts for the majority of stars 
{Ak < 0.02 mag), and thus we have ignored the effects of 
circumstellar extinction. 



4 CCD OBSERVATIONS 
4.1 Equipment 

An automated photometric telescope was constructed from 
relatively low-cost commercial components. A Nikon 80- 
200-mm //2.8 ED telephoto lens {D=77 mm) coupled to an 
SBIG ST8-XE CCD detector with 9 fim pixels was attached 
to a Losmandy G-11 German Equatorial Mount controlled 
by an Astrometric Instruments SkyWalker II stepper motor 
system (Figure |4|. 

The telescope was housed in a roll-off roof observa- 
tory in a suburban backyard in Canberra, Australia, at lati- 
tude 35° S. Observations were gathered on nearly every clear 
night over a 5.5-year period using a fully automated software 
pipeline that controlled target selection and scheduling, tele- 
scope slewing, closed-loop pointing corrections, and camera 
operations (V. Tabur, PhD thesis, in prep.). In addition to 
the photometric survey described here, the observatory has 



4.2 Detector Calibration 

The detector bias, readout noise and gain were determined 
experimentally and found to be in close agreement with the 
manufacturer's stated values. 

We used the procedure outlined by |Berry fc Burneli] 
( |2000 1 to determine that our detector was linear below 47000 
ADU to within 0.1%. Stars with pixel values exceeding this 
limit (prior to calibration) were marked as saturated and 
were discarded by our photometric pipeline. Integrations 
shorter than 2 s were found to be non-linear due to a 0.02 
s delay in closing the mechanical vane shutter. A shutter 
mapping correction was derived using the method of IZissell] 
(2000) but was found to be unnecessary in practice because 
integrations longer than 3 s were used, to reduce scintillation 
noise. 



4.3 Reducing Flux 

Three methods were employed to avoid detector satura- 
tion. By stopping-down the effective aperture of the lens 
to 50 D ^ 64 mm (during different phases of the survey) 
we reduced the number of photons reaching the detector. 
Secondly, short integration times of 3 s ^ tint ^ 25 s were 
employed, with the minimum practical limit set by scin- 
tillation noise and shutter effects. Broad-band photometric 
filters were not used because the SBIG CFWIO filter wheel 
does not provide sufficient back-focus when used in conjunc- 
tion with a camera- lens adaptor. The unfiltered detector re- 
sponse is shown in Figure|5]and, for comparison, the relative 
transmission of a Johnson V^-band filter, and the spectral en- 
ergy distribution of a star of type M3III ( Pickles|1998 1 . 

The third way to reduce saturation was to distribute 
the fiux over a large number of pixels. Stellar images were 
broadened by defocusing the optics, which also helped to 
improve precision by reducing sensitivity to the gate struc- 
ture of the front-illuminated CCD. However, defocusing pro- 
duced a complex PSF due to the multi-element optical de- 
sign, and caused some blending in crowded fields, making 
it difficult to find isolated comparison stars and requiring 
more complicated analysis. 



4.4 Automation 

Once initialised each night, the telescope operated com- 
pletely autonomously, allowing a very efficient survey. The 
observing cycle consisted of the following tasks: 

(i) Target selection: In general, fields were ordered by in- 
creasing R.A. so that those about to set were observed first. 
Allowances were made for exclusion zones to prevent fre- 
quent meridian flips when targets were in the overlap zone 
(visible from both sides of the mount), and to avoid obsta- 
cles. Short slews were weighted more highly to increase the 
duty-cycle. 



(ii) Open-loop slewing: A telescope pointing-model ( Wal- 
lace fc Tritton 1979 1 was used to compensate for mechan- 



ical imperfections such as polar axis misalignment, non- 
orthogonality of the axes, and declination axis fiexure. Point- 
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Table 1. Survey Phases. Observations were concluded on 2008 
Nov 5. 
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Figure 5. Quantum efBciency of the KAF1602 CCD sensor, rela- 
tive transmission of a Johnson V-band filter (used for PEP), and 
the relative flux of a star of spectral type M3III, as a function of 
wavelength. 



ing errors were < 1 arcmin when slewing to the same side 
of the mount, and < 3 arcmin all-sky. Slews were always 
completed with a final movement toward the W and S to 
remove backlash from the gear train. 

(iii) Closed-loop pointing correction: Residual pointing 
errors were corrected using a micro-slew, with the neces- 
sary ofi'sets determined from a 2-s, low-resolution integra- 
tion. The plate solution was calculated using the Optimistic 



Pattern Matching algorithm ( Tabur 2007 1 , which was specif- 



ically designed to handle wide fields of view containing thou- 
sands of stars. The additional integration helped to flush 
the detector of any potential Residual Bulk Image caused 
by near saturation of the detector ( Rest et al.||2002" l, while 
the final pointing correction reduced the effects of residual 
flat-fielding errors ( [Everett fc Howell||1975 1. 



4.5 Survey Phases 

Our survey consisted of three phases (Table [T]). Initially, a 
focal-length of / = 142 mm was employed with integrations 
ranging from 1 s ^ iint ^ 25 s, and with the image slightly 
defocused. The detector was operated at full resolution with 
sub-frames read out from the central 1/3 of the CCD. The 
brightest targets required short exposures of ~ 1-2 s, which 
led to a high level of scintillation noise. Despite the short 
integrations, some data were saturated and had to be dis- 
carded. 

Phases 2 and 3 consisted of greater levels of defocus, 
plus an increased focal length to offset the blending of stellar 
sources. The sub-frame size was enlarged to increase the 
probability of finding bright comparison stars, which were 
required for the ensemble photometry. 

Phase 3 used strongly defocused images yielding Full- 
Width at Half-Maximum (FWHM) values up to 25 pixels. 
Despite the difficulties this caused during analysis, requiring 
blended stellar sources to be deconvolved, this approach was 
adopted because it permitted the use of 25-s integrations for 
all but the very brightest targets, thus reducing scintillation 
noise to ~ 0.5%. Although this could have been reduced 
further with longer integrations, the limit was set to allow 
~125 fields (one visible hemisphere) to be visited within a 
5-hour observing window. 

Multiple consecutive integrations were obtained for each 
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target field to measure internal precision and reduce scintil- 
lation noise. A series of 3, 5, 7, or 9 images were acquired for 
each field utilising 25, 12, 6, or 3 s integrations, respectively. 
The images were analysed separately and instrumental mag- 
nitudes combined to yield mean nightly values for all targets 
and their comparison stars. 

Even the Phase 3 regime was insufficient to prevent the 
two brightest stars, 7 Cru and 13 Gru, from saturating the 
detector when observed near the meridian, although some 
useful images were obtained at higher airmass. Their light 
curves were supplemented with PEP data (Section [6|. 



5 CCD PHOTOMETRY 
5.1 Flat-Fielding 

Reliable flat-field images are notoriously difficult to produce 
for wide flelds. Twilight flats have been shown to contain 
gradients of 1% per degree when obtained too far from 
the null-point, where sky- illumination gradients are min- 
imised (Chromey & Hasselbacher 19961, while blank sky 



flats are non-uniform over 0.5° scales due to natural altitude- 
dependent gradients ( Zhou et al.||2004 l. 

We used a combination of approaches to solve these 
problems. Firstly, high S/N dome flats were acquired 
monthly to model the high-frequency, inter-pixel variations 
in detector response. Illumination gradients of ~ 1% were 
present, despite placing a two-stage diffusing screen near the 
entrance pupil of the optics ( [Zhou et al. 2004 1. To remove the 
gradients, each dome flat was corrected to match the low- 
frequency variations present in a master sky-flat (created by 
median combining 59 images obtained on a moonless night 
with the telescope pointing at the zenith) , modelled using a 
3rd-order, 2-dimensional polynomial. Finally, since the uni- 
formity of the master sky-flat was not assured, we used the 
fact that an imperfect flat-field calibration will lead to spa- 
tially dependent photometric errors to determine an illumi- 
nation correction. Since we observed each fleld hundreds of 
times, on both the east and west sides of the mount, the 
180° rotation of the detector that occurs after a meridian- 
flip was used to model the change in mean magnitude of each 
comparison star as a function of its spatial location on the 
detector, analogous to the dithering of stars across an im- 



age to build a Photometric Super Flat (see Manfroid||1995 
Boyle et al.|2003[ ). In this way, we measured gradients in the 



master sky-flat of up to ± 6 mmag, which we corrected in 
the photometric pipeline. 
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5.2 Atmospheric extinction 

To a second-order approximation, the true magnitude of a 
star, Vo, measured outside of the Earth's atmosphere, is re- 
lated to the instrumental magnitude, Finst, by 

Vo = Vinst - k'X - k"XC, (1) 

where X is the airmass, C is the colour-index of the star and 
k' and k" are the first- and second-order extinction coeffi- 
cients. By using long-term observations of the comparison 
stars gathered over a wide range of airmasses, we deter- 
mined mean values for k' and k" to be 0.16 mag airmass"^ 
and 0.024 mag per airmass per unit colour, respectively. We 
applied the correction in the pipeline. KV — I colour index 
was used for the colour term, to match the wide spectral 
response of the unfiltered detector, which spans 450 nm at 
FWHM and extends to the / band where M giants radiate 
much of their energy. 

5.3 Photometric precision 

Magnitudes of targets and comparison stars were deter- 
mined using aperture photometry. Faint, defocused stars 
in the sky-annulus of some comparison stars, particularly 
at low galactic latitudes, caused a systematic increase in 
sky-background measurements during bright moonlit peri- 
ods, when they became difficult to separate from the real 
sky-background using a sigma-clipping routine. The prob- 
lem was solved by mapping the positions of these stars and 
excluding their pixels from the sky annulus. 

The photometric precision, measured as the rms scatter, 
is shown in Figure [6] for 4979 comparison stars. Our target 
stars have instrumental magnitudes brighter than 9, so they 
all fall in the more precise tail of the distribution. The best 
photometric precision achieved is limited to ~ 6 mmag, even 
for the very brightest stars, which indicates the scintillation 
noise due to the small aperture and short integration times. 



Young ( 1967 1 described the well-known relation for atmo- 



spheric scintillation, S, the ratio of rms intensity fluctuation 
to mean intensity, as 



0.09d"^''^X^''^ exp(-/l/8000)/^ 



(2) 



where d is the telescope aperture in cm, X is the airmass, 
h is altitude of the observing site in metres, and tint is the 
integration time in seconds. Substituting values of d = 6.4 
cm, h — 641 m, and tint = 3 x 25 s for observations con- 
ducted during Phase 3, and assuming a mean airmass of 
X=1.5, we estimate the contribution due to scintillation to 
be ~ 7 mmag, which agrees well with the lower limit on the 
measured precision. 

5.4 Ensemble photometry 

Up to 50 of the brightest comparison stars in each field were 
combined to form an ensemble. We assumed constancy of 
the ensemble average, that is. 



(3) 



J2 "nuwi 

where rrii are the instrumental magnitudes of the compar- 
ison stars, Wi are their weights, and ZP is the zero-point 
of each image in the time-series. The ensemble was used to 
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Figure 6. Individual comparison star scatter after correction of 
sky background estimates in crowded fields. Each measurement 
is the mean of three 25 s integrations with airmass < 2.0 (4979 
stars). 



determine iteratively the ZP for each image ( |Balona|1995l ), 
permitting the scatter of each comparison star to be de- 
termined, unstable stars to be excluded, and weights to be 
assigned to the remainder. Following [Gilliland et al. (19931, 
weights were assigned using the reciprocal of their variance, 
normalised to a mean weight of 1, and clipped to a maxi- 
mum value of 2 to prevent bright stars with low-amplitude 
variability from perturbing the solution. 

Aperture photometry in crowded fields was compro- 
mised by the strong defocus, which caused some blending of 
stellar sources. Comparison stars were generally unaffected, 
since they were manually selected to avoid nearby compan- 
ions. In cases where some contamination was unavoidable, 
the stars exhibiting the largest errors were excluded from 
the ensemble. 

Figure [T] shows the distribution of the uncertainties for 
our ensemble of comparison stars, which were mostly about 
3-5 mmag. Since the M giants were typically the brightest 
stars in each field, their internal errors were small, result- 
ing in a total uncertainty for their differential magnitude 
determination of 7-8 mmag. 



5.5 Optimal Image Subtraction 

Nine M giant target^ were blended with nearby stars and 
required a Difference Image Analysis technique to subtract 
the fiux contributed by their companions, allowing the vari- 
able sources to be measured free of contamination. Following 



Alard & Lupton ( 1998j), we implemented an Optimal Image 
Subtraction (OIS) pipeline. High S/N images taken under 
the best seeing conditions were selected as reference images. 
Due to a small, spatially dependent asymmetry in the de- 
focused PSF, two reference frames were selected for each 
field, corresponding to pre- and post-meridian-flip orienta- 
tions. A high-order 2D polynomial was used to register each 
image with its reference to cater for the presence of optical 



1 HIP 40144, 44281, 46027, 64117, 73764, 80047, 83258, 87390, 
101810 
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Figure 7. Distribution of mean ensemble uncertainty for 243 
survey fields using 1 milli-mag bins. 



distortions caused by the fast (//2.8) optics. Well-focused 
images yielded residuals of ~0.05 pix, but the Phase 3 im- 
ages had larger residuals (~0.15 pix), due to the difficulty 
in determining accurate centroids for the defocused PSFs. 
Since our survey fields were not crowded, we followed 



Kochanski et al. ( 1996 1 by only using the pixels associated 



with the 50 brightest stars, to prevent the Poisson noise of 
the blank sky background from contributing to the kernel 
solution. Since the tuning of OIS parameters is not straight- 
forward, we adopted the method and recommendations of 



Israel et al. (20071 to determine optimal values. 



We found that our OIS pipeline greatly improved the 
photometry in crowded fields. The noise in the power spec- 
tra of hght curves produced using OIS were reduced, as were 
the strengths of spurious peaks at periods of 1 and 1/2 yr 
(Figure [8|. We used two Cepheids (AT Pup and RS Pup), 
located in the field of V427 Pup, as test cases. Both have 
well-known periods and AT Pup is strongly blended with 
a brighter companion. Figure |9] shows their phased light 
curves, and demonstrates the dramatic improvement in pho- 
tometric precision when using OIS on blended objects (up- 
per panel). In some cases OIS did not produce optimal re- 
sults, such as the slight increase in scatter for RS Pup, which 
we attribute to sampling issues and numerical instabilities. 
We restricted its use to 9 blended targets. 



6 PHOTOELECTRIC PHOTOMETRY 

Independent PEP data for 23 of the project stars were ob- 
tained by one of us (T.T.M.) as part of a survey of south- 
ern SRVs ( |Moon, Otero fc Kiss|[2007| . The PEP observa- 
tions span a similar temporal range to the CCD survey, but 
are generally sparser due to the manual operation of the 
equipment. Each measurement is the mean of 5 consecu- 
tive 10-second integrations using V and B broadband filters, 
obtained in the sequence Vstar, Bstar, Bsky, V^aky, and brack- 
eted by observations of two comparison stars. Measurements 
were corrected for extinction and transformed to standard 
V magnitudes and B — V colour indices using data from the 
General Catalogue of Photometric Data ( [Mermilliod et aT] 
T997|. 



6.1 Transformation to Vj 

In order to combine CCD and PEP data, the unfiltered CCD 
observations of the 23 stars were transformed to Johnson V 
magnitudes using numerous comparison stars in each field. 
Since some comparison stars were likely to be variable, en- 
sembles were used to identify stable stars and their mean 
instrumental magnitudes. In the absence of colour informa- 
tion in the instrumental system, the {V — I)c colour index 
of each comparison star was used to fit the relation 



Vj - Vins = Co./ +Cl{V^I)+ C2{V - If 



(4) 



Here, Vj is the Johnson V-band magnitude, Vins is the mean 
instrumental magnitude, cq,/ is a zero point for each field, 
and ci and C2 are the colour terms. Separate zero-points 
were required for each field because the ensemble used to 
derive the mean instrumental magnitudes consisted of an 
inhomogeneous set of stars and weights. 

After fitting, some stars had an excessively large differ- 
ence between the catalogued Vj and the transformed value, 
caused by contamination of their photometric aperture by a 
nearby companion, resulting in the asymmetry in Figure [Tol 
Outliers were removed using 3 iterations of a sigma clip- 
ping routine where stars > 2.5a from the median residual 
were removed. 4562 stars were retained in the final list of 
well-fitting, uncontaminated comparison stars, resulting in 
a scatter of 0.06 mag. 

Although the fit was reasonable for the comparison 
stars, the cooler M giants exhibited larger residuals with 



systematic differences between the two systems (Figure 11 
left panel). We attribute these to the lack of instrumental 
colour information for use in the transformation, together 
with the use of an unfiltered detector that was sensitive to a 
broad range of wavelengths, and extrapolating the relation 
to cool temperatures beyond the range of comparison star 
colour indices {V ~ I > 1.7). The left panel shows an ap- 
proximately linear systematic difference between Vj in the 
CCD and PEP systems for the M giants, as a function of 
spectral type. The filtered PEP data are approximately 0.05 
mag brighter than the unfiltered CCD for M2III stars, and 
the difference increases to 0.37 mag for later- types. A hint 
of the same trend is visible in Figure [To] for 2 < V — I < 2.5, 
corresponding to spectral types up to M2III, though there 
are too few stars to be certain, since most M giant compar- 
ison stars were variable and were thus excluded from the 
fit. 

To adjust for the difference, a zero-point correction was 
added to align the CCD and PEP fight curves. The correc- 
tion was determined by using the CCD observations that 
bracketed each PEP observation by ±2 days to interpolate 
the CCD magnitude to the same date as the PEP obser- 
vation. The mean difference between the CCD and PEP 
magnitudes defined the offset. 

(right panel) shows the rms uncertainty of 
V^j.PEp) residuals after applying the zero-point 



Figure 

(Vj-,CCD 

correction, as a function of spectral type. Most stars have 
a small scatter of ~ 0.02 - 0.04 mag, somewhat larger than 
the uncertainty of the PEP measurements, which were gen- 
erally better than 0.01 mag. We attribute the differences 
to errors in the transformation. However, as expected, the 
coolest stars exhibit a larger scatter, caused by a systematic 
increase in passband difference due to the transmission of 
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Figure 8. Light curve and power spectrum of NSV 20880 (HIP 83258), which was blended with a bright companion, before (top) and 
after processing with OIS (bottom), resulting in the removal of the 1- and 0.5 yr peaks which dominate the original power spectrum. 
Note the change in scale of the corrected power spectrum. 




Figure 9. Phased light curves for the Cepheids AT Pup (top) 
and RS Pup (bottom), showing raw differential magnitudes (red 
pluses) and those following the application of OIS (blue crosses). 
AT Pup is blended with a brighter companion, leading to sig- 
nificant seeing-dependent aperture contamination, while RS Pup 
does not have a bright companion. 



longer wavelengths on the red side of the l^-band cut-ofT, to 
which our unfiltered CCD is sensitive. 

The PEP and CCD hght curves exhibit the same gen- 
eral features, although the former exhibit greater V-band 



amplitudes (Figure 12 1. The most extreme example occurs 
for the M8III variable V744 Cen, which has a 36% greater 
amplitude in the PEP system. Fortunately, the coolest stars 
exhibit long-period, large-amplitude variations and so their 
periods are well determined using CCD data alone (Table 
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Figure 10. The difference between catalogued Vj and the trans- 
formed values as a function oi V — I colour index. Stars with 
an excess are involved with a close companion that increased the 
measured flux within the photometric aperture, resulting in the 
asymmetrical distribution. 



Table 2. Late- type stars with large residuals after fitting a zero 
point to transformed CCD observations. PEP data were ignored 
when determining periods for these stars. 



HIP 


Name 


Amp 


Std Dev 


Period 


Spectral 






(mag) 


(mag) 


(days) 


Type 


61404 


BO Mus 


0.55 


0.069 


133 


M6II/HI 


66666 


V744 Cen 


1.11 


0.122 


167 


M8III 


98608 


u Pav 


0.45 


0.066 


87.8 


M6III 


106044 


SX Pav 


0.54 


0.064 


50.5 


M5III 



6.2 VZ Cam 

VZ Cam was observed by one of us (B.Sz.) using the 24" 
telescope of the Konkoly Observatory between 1968 and 
1977. The data were obtained with an RCA 1P21 photo- 
multiplier through Johnson UBV filters, using two compar- 
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II Cepheid RU Cam ( Szeidl et al. 1992 1 . We added the object 
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Figure 11. Mean difference between M-giant Vj magnitudes 
in the CCD and PEP systems as a function of Morgan-Keenan 
spectral type (left), and the resulting Icr scatter in the differences 
after applying a zero-point correction as described in the text 
(right). Late- type stars exhibit greater errors due to increasing 
differences in effective passband. 



ison stars (comp^HD 57508, 1^=6.550, B - V^Q.mQ and 
U - 5=0.63; check=HD 55075, V'=6.325, B ~ l/=-0.040 
and U — B—~0.05). The typical photometric accuracy is 
±0.01-0.02 mag. The star was observed in a program for 
obtaining long-term photometric data for the peculiar Pop. 



to our sample because of its accurate Hipparcos parallax and 
the previously unknown two periods that can be determined 
from the decade-long light curve. This is the first analysis of 
VZ Cam in the literature. 



7 PERIOD ANALYSIS 

Characteristic periods were determined for each fight curve 
using a weighted Discrete Fourier Transform (DFT). The 
weighted mean magnitude of the set was first subtracted to 
prevent a peak at zero frequency from dominating the spec- 
trum. Observations were weighted by the inverse square of 
their uncertainty, except those combined with PEP data, 
which were analysed using multiple weighting schemes as de- 
scribed below. The DFT code was an implementation of the 



algorithm proposed by Scargle (1989 Appendix D), modi- 
fied to include weights for each point. 

Observations were predominantly spaced at intervals of 
1 day, resulting in a Nyquist frequency of 0.5 cycle day~^. 
Power spectra were created for the frequency range 0-0.4 
cycle day~^. When searching for the most significant peak 
(greatest power), only frequencies < 0.2 cycle day~^ were 
considered, corresponding to a pulsation period > 5 d. The 
mean power in the frequency range 0.2-0.4 cycle day~^, 
where no pulsation signal is expected, was used to estimate 
the noise level. The peaks at low frequency have lower S/N 
because of the rising noise fioor, caused by 1// noise. Only 
periods shorter than 300 d were kept, since we consider them 
to be the most reliable. Longer periods were less well de- 
termined due to gaps in the sampling, possible instrumental 
drifts, and difficulty matching light curves in different phases 



due to insufficient comparison stars in the earlier, smaller 
FOVs. 

Light curves containing both CCD and PEP data 
showed higher noise in their corresponding DFTs compared 
to those prepared using CCD data only. Noise levels were at 
a minimum when all observations were weighted equally, but 
increased for weighting schemes using and a~^. This ef- 
fect was caused by the l/-filtered PEP observations having 



larger amplitudes than the unfiltered CCD data (Figure 12 1, 



probably caused by the lack of instrumental colour informa- 
tion in the transformation. However, since the periods are 
well-correlated, we find that despite the increase in noise, 
the fitted periods remain unaffected. Since equal numbers 
of PEP and CCD observations tend to produce the greatest 
noise, the weighting scheme was selected dynamically based 
on the relative number of observations. An inverse variance 
was used when the number of CCD observations exceeded 
the PEP observations by a factor of two, otherwise equal 
weights were used. 

Periods were extracted using iterative sinewave fitting 
( Frandsen et al.||1995 l. That is, least-squares minimisation 
was used to fit a sinusoid with a frequency corresponding 
to the strongest peak in the DFT, while solving for the op- 
timum amplitude and phase. The best-fitting sinusoid was 
subtracted from the data, and the process repeated until no 
strong peaks (S/N 7.5) remained. 



7.1 The catalogue 

The survey result^ are presented in Table [sj with stars or- 
dered by ascending Hipparcos catalogue number (column 1). 
Columns 2 and 3 list the star name, and spectral type from 
the Bright Star Catalogue ( jHoflleit fc Warren||T99l| ) where 
available, otherwise from the General Catalogue of Variable 
Stars ( Samus et al.|[2009 GCVS), or the most recent entry 
from the compilation of spectral types by Skiff ( 2009| ) . Addi- 
tional columns list the Johnson l^-band magnitude, revised 



Hipparcos parallax and uncertainty (van Leeuwen 2007[ ) 



apparent K magnitude, absolute K magnitude and uncer- 
tainty, if-band magnitude source (described in Section [3] 
D=DIRBE, G=Gezari, 2M=2MASS), number of CCD and 
PEP observations, and detected periods. 

GZ Eri and NSV 16260 are fisted as G and K giants, 
respectively, although other sources list them as M giants 
( jHouk fc Smith-Moore|l9 88 Houkj|l982|. The shape of their 
light curves and the derived periods are consistent with the 
usual variability in early M-type stars. Fourteen stars have 
no determined periods, although we retain them in Table [S] 
to document their Mk- Five of these stars (BW Oct, NSV 
18694, V931 Her, V983 Her, and V2158 Cyg) have no obser- 
vations due to physical obscuration or low altitude. A further 
6 stars (7 Hyi, A Vel, V1919 Cyg, V2090 Cyg, V2114 Cyg, 
V345 Peg) have insufficient observations for reliable period 
determination. The remaining 3 stars have either long peri- 
ods, or periods with low S/N, both of which were excluded 
by our selection criteria, as follows. X Tra shows two periods 
of ~ 385 and ~ 455 d. CI Boo shows a long period of ~ 485 
d and a shorter period of ~ 51 d, although the latter is only 
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Figure 12. CCD (top panel) and PEP light curves of f3 Gru. The agreement is good, with the PEP observations exhibiting a slightly 
larger photometric amplitude. 



detected at S/N ^ 7.2. V440 Aur has a possible 39 d period 
which is only detected at S/N ~ 5.8. 



7.2 Sample light curves 

A selection of light curves and their corresponding power 
spectra are shown in Appendix A in the online version of 



the paper (see Supporting Information). Figures Al and A2 



show an example light curve and power spectrum for Z Eri. 
The full set of plots are available on request from the au- 
thors, and photometric data are available electronically from 
CDS, Strasbourg. 

Relatively few stars show Mira-like light curves with 
regular, high-amplitude, mono-periodic behaviour, although 
our selection in the appendix is biased toward higher- 
amplitude variability for aesthetic reasons (see for example, 
DU Cha and AY Dor). Instead, the majority of stars exhibit 
short-period variability on time-scales of 10-100 d, small 
amplitudes {A < 0.2 mag), and multiple periods, resulting 
in complicated light curves (see, for example, V512 Car). 

It is well known that the excited frequencies in SRVs are 



not stable over time, and may vary by a few percent | |Mattei 
et al.||1997[ [Kiss et al.||1999| ). The amplitudes of the modes 



can vary on timescales of 1-5 yr (Kiss et al. 2000 Percy 



|et al.|p003| , and in some cases, phase variations have been 
observed on short time-scales, consistent with the properties 



of a stochastically excited and damped oscillator ( Bedding et 
|al.|2005| . Consequently, the complex light variations of SRVs 
cannot be fully modelled by a simple sum of sine waves. The 
occurrence of closely separated frequencies is quite common 
amongst our sample (see, for example, Z Eri, WX Hor, e 
Mus, BL Cru), and is probably caused by modulations in 
the amplitude and/or phase of a single oscillation mode, and 
may indicate the finite lifetime of the oscillation. We did not 
try to combine periods by introducing arbitrary criteria, but 
instead, provide the raw data in Table |3] 

We found that 63% of periods are shorter than 50 d and 
nearly 5% of periods are shorter than 15 d. Of these, the 
shortest reliable period is 8.7 d (NSV 15560, A —9 mmag) 



and corresponds to the strongest peak in its power spectrum 
(see Fig A6). The shortest period belongs to NSV 18346 
(P=6.7 d, M3III, J - K=1.16, A'/k=-4.7) but has an am- 
plitude of only 5 mmag and may be spurious. 

Although we did not retain periods longer than 300 d 
m Table m we note that WW Pic exhibits a long period of ~ 
370 d and a short period of 36.1 d. The period ratio of ~ 10 is 
consistent with galactic SRVs having long secondary periods 
(LSP). Although the precise nature of their variations are 
still unknown ( jOliver fc Wood|2003| [Wood et al.|200"4l ), LSP 
behaviour has been identified amongst SRVs in the lower 



metallicity environment of the LMC (Wood et al. 19991. 
Moreover, WW Pic has an Mi<:=-4.74, placing the 370- 
d period on sequence D of [Wood et al.| ( |1999[ ) , which is 
composed of LMC stars having LSPs. 
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11 


(0 


38) 


D 


405 





97.9 116.4 193.4 225.7 


57613 


H Hya 


M4+III 


5.10 


5 


98 


(0 


27) 


-0 


30 


— 6 


43 


(0 


10) 


D 


381 





29.5 30.4 39.5 41.4 89.1 


57800 


RU Crt 


M3 


8.14 


7 


57 


(0 


98) 


+0 


07 


— 5 


54 


(0 


28) 


D 


328 





60.6 66.4 101.1 


59588 


V335 Hya 


M4III 


6.34 


2 


79 


(0 


46) 


-0 


17 


-7 


96 


(0 


36) 


D 


429 





107.6 141.2 


59929 


e Mus 


M5III 


4.06 


10 


82 


(0 


17) 


-1 


43 


—6 


26 


(0 


03) 


G 


408 


251 


32.1 32.7 42.5 43.7 44.9 46.0 63.4 196.5 


60333 


KG Vir 


M6 


7.64 


4 


53 


(0 


65) 


+1 


71 


—5 


02 


(0 


31) 


D 


314 





18.4 22.3 24.6 32.1 56.1 116.7 157.0 


60421 


TT Crv 


M3III 


6.48 


3 


53 


(0 


39) 


+1 


64 


— 5 


64 


(0 


24) 


D 


332 





15.4 15.6 19.5 19.9 26.8 37.7 182.8 


60781 


BL Cru 


M4-5III 


5.38 


7 


36 


((3 


32) 


-0 


19 


— 5 


87 


(0 


10) 


D 


482 


179 


3(3.7 42.3 43.6 


60979 


V928 Ccn 


M2II-III 


6.(3(3 


4 


42 


((3 


37) 


+1 


16 


— 5 


63 


(I) 


19) 


D 


434 





19.(3 19.2 22.3 24.6 29.6 46.4 102.2 


61084 


7 Cru 


M3.5III 


1.59 


36 


83 


((3 


18) 


-3 


10 


— 5 


27 


{(> 


06) 


G 


209 


232 


12.1 15.1 16.5 54.8 82.7 104.9 


61404 


BO Mus 


M6II-III 


6.11 


2 


97 


(0 


32) 


-0 


83 


— 8 


48 


(0 


23) 


D 


488 


156 


134.8 229.9 


61658 


FW Vir 


M3+III 


5.68 


6 


62 


(0 


31) 


+0 


94 


—4 


96 


(0 


10) 


D 


301 





9.3 28.6 89.7 


61908 


V341 Hya 


M3III 


6.64 


4 


32 


(0 


43) 


+1 


72 


—5 


12 


(0 


23) 


D 


415 





8.6 18.3 18.7 36.2 216.6 281.7 




LM Mus 




7.11 




87 


(0 


47) 


+1 


40 


— 5 


18 


(0 


32) 









2U.cS ^o.y OU.U OU.O 


62985 


!/> Vir 


M3-III 


4.77 


5 


99 


(0 


23) 


+0 


14 


-5 


99 


(0 


09) 


D 


348 





22.4 23.5 24.5 30.1 31.3 49.4 162.6 


63090 


S Vir 


M3+III 


3.39 


16 


44 


(0 


22) 


-1 


24 


-5 


16 


(0 


03) 


D 


229 





13.0 17.2 25.6 110.1 125.8 


63752 


KX Vir 


M4.5III 


7.41 


3 


19 


((3 


71) 


+1 


76 


— 5 


73 


{(> 


49) 


D 


332 


(3 


15.6 17.8 23.4 25.6 69.4 77.6 250.0 


63950 


FS Com 


M5III 


5.53 


4 


43 


((3 


41) 


-0 


20 


-6 


97 


{(> 


20) 


D 


253 





33.7 44.6 46.8 159.5 


64117 


NSV 06103 


Mini 


6.34 


3 


55 


(0 


44) 


+1 


81 


— 5 


46 


(0 


41) 


2M 


313 





14.4 16.8 25.4 183.8 250.6 


64645 


LO Vir 


M6 


7.30 


3 


75 


(0 


53) 


+1 


64 


-5 


50 


(0 


31) 


D 


332 





20.0 22.0 22.3 26.8 91.7 


64852 


a Vir 


Mini 


4.78 


4 


83 


(0 


19) 


+0 


49 


-6 


10 


(" 


09) 


D 


328 





23.4 24.3 27.9 34.1 36.8 64.0 166.3 


65311 


NSV 06212 


Mini 


6.17 


3 


99 


(0 


34) 


+1 


70 


-5 


31 


(0 


19) 


D 


440 





17.9 21.7 22.4 22.7 31.7 34.0 52.2 99.7 


66666 


V744 Cen 


M6III 


5.74 


6 


35 


(0 


33) 


-0 


72 


-6 


72 


(0 


12) 


D 


461 


75 


96.8 102.2 166.9 


66783 


LY Mus 


M4III 


6.91 


3 


48 


(0 


41) 


+0 


95 


-6 


37 


(0 


26) 


D 


405 





53.1 55.3 84.7 106.6 128.7 261.8 


67288 


NSV 06442 


M2III 


5.41 


4 


82 


(0 


39) 


+1 


01 


-5 


60 


(0 


18) 


D 


376 





18.3 18.5 18.8 25.4 26.9 50.2 


67457 


V806 Cen 


M4.6III 


4.19 


17 


82 


(0 


21) 


-1 


67 


-5 


41 


(0 


03) 


D 


354 


107 


26.6 
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HIP 


Name 


Spectral 
Type 


my 
(mag) 






7r 

(mas) 


rriK 
(mag) 






Mk 
(mag) 


Src 


Nc 


Np 


Periods 
(days) 


68815 


6 Aps 


M7III 


6.69 


8 


84 


(0 


49) 


-2 


02 


—7 


29 


(0 


12) 


D 





259 


100.9 111.0 


68937 


ER Vir 


M4III 


6.57 


4 


16 


(0 


50) 


+ 1 


45 


— 5 


48 


(0 


26) 


D 


372 





16.1 18.8 25.7 29.6 80.2 


69269 


ET Vir 


M2III 


4.93 


7 


08 


(0 


32) 


+0 


58 


— 5 


18 


(0 


18) 


D 


367 





22.6 23.8 36.4 37.6 39.8 48.8 259.1 


69614 


FS Vir 


M4III 


6.41 


4 


04 


(0 


52) 


+ 1 


54 


— 5 


44 


(0 


28) 


D 


351 





13.5 16.1 17.0 22.6 144.7 170.1 


69850 


NQ Aps 


M3III 


7.28 


2 


82 


(0 


62) 


+ 1 


93 


— 5 


84 


(0 


49) 


D 


241 





21.2 28.2 


70236 


CI Boo 


M3III 


6.50 


4 


23 


(0 


47) 


+0 


58 


—6 


29 


(0 


24) 


D 


141 







72122 


HD 129902 


Mini 


6.06 


3 


46 


(0 


51) 


+1 


78 


—5 


54 


(0 


32) 


D 


283 





11.3 19.8 32.3 211.0 


72432 


V768 Cen 


M3III 


5.89 


6 


34 


(0 


34) 


-0 


44 


—6 


44 


(0 


12) 


D 


460 





41.4 57.4 60.9 145.1 191.2 


72655 


NSV 20194 


MOIII 


7.26 


3 


83 


(0 


67) 


+2 


37 


—4 


73 


(0 


38) 


G 


515 





16.8 19.2 19.4 28.3 47.7 192.3 


73589 


EV Boo 


M4III 


6.60 


3 


08 


(0 


41) 


+0 


70 


—6 


86 


(0 


29) 


G 


137 





31.2 35.3 59.3 


73764 


GM Lup 


M6III 


6.29 


5 


40 


(0 


56) 


+0 


00 


— 6 


35 


(0 


23) 


G 


387 





29.9 43.6 45.5 61.3 80.9 243.3 


74582 


X TrA 


C5,5(Nb) 


5.75 


2 


78 


(0 


44) 


-0 


62 


— 8 


39 


(0 


34) 


D 





160 




74982 


FZ Lib 


M4III 


6.95 


2 


84 


(1 


03) 


+1 


00 


— 6 


76 


(0 


79) 


D 


357 





22.2 25.9 26.4 62.3 


74999 


OO Aps 


M2III 


f).59 


2 


82 


(0 


45) 


+1 


81 


— 5 


93 


(0 


36) 


D 


67 


103 


24.9 


76423 


Ser 


M5IIb-III 


6.51 


4 


86 


(0 


46) 


— 1 


02 


— 7 


60 


(0 


21) 


D 


329 





86.7 173.0 244.5 


76573 


LY Ser 


M4III 


6.82 


4 


75 


(0 


47) 


-0 


46 


—7 


09 


(0 


22) 


D 


243 





66.7 76.5 220.8 


79072 


FS Ser 


M3.5III 


5.69 


4 


04 


(0 


43) 


+0 


56 


—6 


41 


(0 


23) 


G 


354 





20.6 25.2 25.6 26.4 48.4 


79086 


FQ Ser 


M4III 


6.47 


4 


92 


(0 


73) 


+0 


54 


—6 


01 


(0 


33) 


G 


354 





31.3 32.1 36.2 41.4 44.0 71.6 144.1 


79349 


LQ Her 


M4.5III 


5.74 


4 


78 


(0 


31) 


+0 


15 


—6 


46 


(0 


14) 


D 


274 





24.6 25.1 28.0 29.9 33.6 


79754 


V368 Nor 


M2III 


5.45 


4 


21 


(0 


43) 


+0 


62 


— 6 


27 


(0 


22) 


G 


396 





18.9 21.4 22.1 24.2 33.6 


80047 


o Aps 


M5III 


4.68 


4 


28 


(0 


16) 


-0 


73 


— 7 


59 


(0 


08) 


Ct 


307 





68.0 94.9 101.7 


80620 


V2105 Oph 


M3-III 


5.24 


5 


65 


(0 


39) 


+0 


48 


— 5 


80 


(0 


15) 


D 


371 





20.5 21.2 21.4 41.6 206.6 


81376 


V840 Ara 


M2/M3III 


6.84 


4 


67 


(0 


57) 


+1 


84 


— 4 


83 


(0 


38) 


2M 


384 





16.3 17.9 24.3 42.6 51.1 76.6 164.7 


81426 


V904 Her 


M2III 


6.96 


2 


87 


(0 


52) 


+2 


19 


-5 


54 


(0 


40) 


D 


224 





22.8 65.5 95.4 257.7 


82028 


HR 6227 


M3III 


5.60 


7 


19 


(0 


42) 


+0 


82 


—4 


91 


(0 


13) 


D 


346 





13.2 14.7 17.1 21.9 38.7 


83258 


NSV 20880 


M4III 


7.39 


3 


05 


(0 


64) 


+1 


74 


—5 


86 


(0 


46) 


G 


346 





13.3 14.6 15.8 18.9 181.2 


83462 


V931 Her 


M4III 


6.19 


4 


35 


(0 


33) 


+0 


70 


—6 


11 


(0 


17) 


G 










84744 


NSV 21509 


Mini 


6.80 


3 


42 


(0 


48) 


+2 


48 


—4 


86 


(0 


35) 


D 


160 





19.3 


84780 


V2113 Oph 


M5III 


6.43 


6 


43 


(0 


52) 


-0 


21 


-6 


18 


(0 


18) 


D 


372 





47.3 49.3 70.5 76.1 267.1 


84833 


V656 Her 


M2III 


5.(31 


6 


80 


(0 


36) 


+0 


89 


— 4 


96 


(0 


13) 


D 


343 





22.3 93.9 176.1 234.7 283.3 


84938 


NSV 21626 


M2.5III 


6.35 


3 


26 


(0 


49) 


+2 


08 


— 5 


37 


(0 


33) 


D 


348 





19.7 23.6 


85435 


V859 Ara 


M2II 


7.07 


3 


90 


(0 


58) 


+ 1 


94 


— 5 


12 


(0 


40) 


2M 


372 





14.6 18.9 24.9 35.2 97.8 126.3 


85760 


NO Aps 


M3III 


5.83 


3 


69 


(0 


28) 


+0 


75 


-6 


44 


(0 


17) 


D 





207 


26.2 26.6 


85934 


V642 Her 


M4III 


6.36 


5 


42 


(0 


52) 


+0 


95 


—5 


40 


(0 


21) 


D 


364 





21.9 25.5 33.6 35.4 61.2 130.6 163.4 229.9 


86153 


V959 Her 


M4III 


6.37 


4 


23 


(0 


44) 


+1 


41 


—5 


46 


(0 


24) 


D 


188 





26.3 35.7 


87390 


HD 162189 


M2III 


5.94 


5 


77 


(0 


66) 


+1 


80 


—4 


40 


(0 


21) 


G 


207 





14.1 28.5 


88563 


V980 Her 


M2III 


6.97 


3 


39 


(0 


63) 


+2 


27 


—5 


10 


(0 


42) 


2M 


358 





19.1 24.6 26.0 31.4 135.3 


88761 


V983 Her 


M... 


7.33 


3 


28 


(0 


47) 


+2 


38 


-5 


05 


(0 


32) 


D 










89172 


V669 Her 


M3III 


4.96 


5 


99 


(0 


22) 


+ 


33 


— 5 


79 


(0 


08) 


D 


152 





22.9 24.0 


89527 


V2392 Oph 


M4III 


6.10 


4 


22 


(0 


50) 


+ 


68 


— 6 


22 


(0 


26) 


D 


379 





17.6 20.0 25.9 32.0 33.2 97.7 


91135 


V4401 Sgr 


M4III 


6.77 


6 


46 


(0 


56) 


+ 


72 


— 5 


24 


(0 


19) 


G 


387 





18.8 23.1 23.8 24.2 33.4 36.9 


91789 


V536 Lyr 


M... 


8.20 


2 


64 


(0 


71) 


+ 2 


11 


— 5 


81 


(0 


59) 


D 


127 





50.8 282.5 


92079 


V4405 Sgr 


M4III 


6.34 


4 


49 


(0 


40) 


+0 


42 


—6 


33 


(0 


19) 


G 


291 





21.3 26.2 32.6 38.8 66.9 


93270 


V387 Vul 


M3.5III 


6.28 


4 


63 


(0 


41) 


+0 


81 


—5 


89 


(0 


21) 


G 


294 





16.8 19.2 


93919 


NSV 24691 


Ml 


7.53 


5 


45 


(0 


72) 


+1 


93 


—4 


40 


(0 


29) 


G 


380 





12.0 15.5 103.2 


94646 


HD 180377 


M2III 


6.46 


3 


88 


(0 


53) 


+1 


79 


—5 


29 


(0 


30) 


G 


331 





14.9 16.5 16.7 34.6 296.7 


94969 


V389 Vul 


M4 


6.95 


4 


27 


(0 


57) 


+0 


39 


~6 


48 


(0 


29) 


G 


307 





36.1 49.4 58.9 82.3 226.8 


96111 


V4419 Sgr 


M4III 


6.58 


4 


59 


(0 


62) 


+1 


03 


— 5 


68 


(0 


30) 


D 


414 





29.3 39.9 57.4 


96310 


V1919 Cyg 


M3III 


6.64 


3 


17 


(0 


39) 


+1 


50 


— 6 


01 


(0 


27) 


G 


15 







96399 


NSV 24813 


M4/5III 


f).71 


4 


63 


(0 


59) 


+1 


21 


— 5 


47 


(0 


28) 


D 


460 





11.9 12.0 13.9 15.3 21.4 29.0 136.1 


97142 


V2090 Cyg 


Mini 


6.39 


3 


36 


(0 


36) 


+0 


74 


— 6 


65 


(0 


24) 


G 


11 







97365 


<5 Sge 


M2II+A0V 


3.68 


5 


49 


(0 


72) 


-0 


78 


—7 


09 


(0 


28) 


D 


286 





32.4 34.0 34.8 64.2 


98438 


VZ Sge 


M4III 


5.33 


3 


20 


(0 


31) 


-0 


10 


—7 


59 


(0 


21) 


G 


319 





36.6 39.2 51.4 65.1 


98500 


V345 Tel 


M4III 


6.75 


5 


40 


(0 


59) 


+1 


50 


—4 


85 


(0 


24) 


G 


461 





17.2 20.5 22.4 42.3 75.4 


98608 


NU Pav 


M6III 


4.95 


6 


86 


(0 


26) 


-1 


47 


— 7 


29 


(0 


09) 


D 


434 


112 


59.0 62.0 87.0 124.7 272.5 


98688 


V3872 Sgr 


M4III 


4.43 


7 


27 


(0 


18) 


-0 


76 


— 6 


46 


(0 


06) 


D 


367 





24.0 30.4 31.3 42.8 60.6 234.7 


98954 


V1472 Aql 


M2.5III 


6.37 


7 


92 


(1 


07) 


+1 


73 


— 3 


78 


(0 


29) 


G 


318 





100.2 


100051 


V2114 Cyg 


M5 


7.23 


2 


50 


(0 


51) 


+0 


86 


— 7 


18 


(0 


44) 


G 


12 







101316 


MT Del 


M4II-III 


7.12 


3 


39 


(0 


66) 


+2 


13 


— 5 


23 


(0 


42) 


G 


307 





26.0 29.0 50.7 64.1 103.7 


101810 


EU Del 


M6III 


6.22 


8 


56 


(0 


50) 


-1 


12 


—6 


46 


(0 


14) 


G 


266 





44.0 60.8 67.3 132.6 235.3 


102096 


AV Mic 


M3II 


6.28 


2 


85 


(0 


40) 


+1 


10 


—6 


64 


(0 


31) 


D 


438 





22.3 23.3 30.3 31.0 32.3 46.0 110.7 


102624 


EN Aqr 


M3III 


4.43 


5 


57 


(0 


28) 


-0 


22 


—6 


50 


(0 


11) 


D 


346 





20.2 24.9 27.2 35.0 36.9 143.9 197.2 


102720 


FI Vul 


M3 


7.65 


5 


15 


(0 


60) 


+0 


16 


— 6 


29 


(0 


25) 


D 


270 





53.8 64.6 73.7 111.7 


103645 


V2142 Cyg 


M3III 


6.85 


3 


53 


(0 


43) 


+2 


47 


— 4 


81 


(0 


30) 


D 


88 





22.3 


104755 


NSV 25504 


M1-2III 


5.06 


3 


65 


(0 


67) 


+ 


69 


— 6 


51 


(0 


40) 


D 


424 





15.5 88.3 130.9 162.6 228.3 


104974 


NSV 13620 


M3III 


5.31 


5 


33 


(0 


33) 


+ 


54 


— 5 


84 


(0 


14) 


D 


348 





19.6 20.1 24.3 91.5 


105334 


T Ind 


C7,2(Na) 


6.15 


1 


72 


(0 


40) 


+ 


60 


— 8 


22 


(0 


51) 


D 





51 


167.6 


105518 


NSV 25556 


M2III 


7.10 


2 


93 


(0 


57) 


+2 


81 


—4 


87 


«o 


49) 


2M 


348 





16.6 


105678 


Y Pav 


C7,3(N0) 


6.28 


2 


48 


(0 


59) 


+0 


33 


—7 


70 


(0 


52) 


D 





152 


232.0 




Viioo v-'ys 


M 


7 77 




18 


(0 


69) 


+1 


99 


4 


92 


(0 


33) 




Q 


Q 




106044 


SX Pav 


M5III 


5.47 


7 


92 


(0 


28) 


-0 


71 


-6 


22 


(0 


34) 


2M 


422 


236 


36.0 50.6 53.1 77.6 88.9 


106062 


NV Peg 


M4.5III 


5.84 


4 


90 


(0 


32) 


-0 


16 


-6 


72 


(0 


14) 


D 


261 





44.0 45.4 77.0 91.5 


107516 


EP Aqr 


M8III 


6.52 


8 


80 


(0 


63) 


-1 


53 


-6 


81 


(0 


16) 


D 


324 





81.8 117.9 125.8 


109212 


OY Peg 


Ml 


6.31 


3 


26 


(0 


38) 


+1 


62 


— 5 


83 


(0 


26) 


D 


260 





28.1 29.4 38.6 41.0 64.3 202.0 


110256 


BO Oct 


M5III 


5.09 


11 


22 


(0 


23) 


-1 


16 


— 5 


92 


(0 


05) 


D 





315 


18.0 23.1 23.4 201.6 


110346 


PT Peg 


M4III 


6.46 


4 


43 


(0 


48) 


+1 


35 


-5 


43 


(0 


24) 


D 


163 





164.7 


110364 


DN Oct 


M3III 


7.30 


2 


89 


(0 
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Figure 13. A comparison of periods for stars appearing in [Glass] 
I&: van Leeuwen] | |2007[ l and this study. Periods are shown for two 
S/N limits. Overlapping symbols indicate that no closer period 
was found by relaxing the S/N limit. 



7.3 Comparison with known periods 

We compared our derived periods to those found in the Ut- 
erature and found that 21 stars in our sample are listed in 



Glass & van Leeuwen (20071. Using each listed period less 



than 300 d, we searched for the closest corresponding period 
from our study. Figure [13] shows a comparison using periods 
selected with two different S/N limits. Although the agree- 
ment is quite good with a S/N limit of 7.5, reducing the 
limit to S /N ^ 5 found matching periods for some lower lu- 
minosity stars, which tend to have smaller amplitudes. The 
agreement is remarkably good considering that SRVs show 
long-term amplitude modulation or mode-switching (Bed- 
ding fc Zijlstra|1998[|Kiss et al.|19"99l ). |Groenewegen| ( [2004 1 
compared the periods of 370 stars over a time span of 17 
years and found that 30 stars showed > 10% change in pe- 
riod, while an additional 36 stars exhibited period changes 
large enough to require assignment to a different P-L se- 
quence. 



7.4 Variability and multi-periodicity 



In their recent review article, Eyer & Mowlavi ( 2008 1 stated 



that "ai/ red giants with spectral types later that early K are 
variable" , which led us to investigate whether we could con- 
firm this within our sample. Excluding 11 stars that lacked 
sufficient observations for period determination, we find that 
all stars are variable at a detection threshold of S/N 5. 
Figure [l4| shows the number of periods detected as a per- 
centage of the effective sample size (250) for a range of S /N 
thresholds. The number of stars with no detected period was 
0, 3, 12, 48 for S/N limits of 5, 7.5 (used for Table|3|, 10, and 
15, respectively. Since we selected the majority of our sample 
using the Hipparcos catalogue's course variability flag (H6; 
section [2|, our results confirm it is a reliable predictor of 
variability. However, several stars exhibited clear variability 
but with amplitudes < 0.01 mag, far less than the ampli- 
tude observed by Hipparcos (H6 2 indicates observed am- 
plitudes ^ 0.06 mag), although it should be noted that the 
amplitudes from iterative sinewave fitting are smaller than 



the peak-to-peak amplitudes of each time-series. We defer 
further discussion of amplitudes to Paper H. 

It is well known that semi-regular variables in the Milky 
Way exhibit multi-periodic behaviour. jMattei et al. (19971 



analysed long-term visual light curves of SRVs using data 
from the AAVSO database, concluding that most SRVs are 
multi-periodic, and exhibit period ratios 1.7 ^ P1/P2 ^ 
1.95. They also confirmed 6 SRVs to be triply-periodic. An 
analysis of long-term visual observations of 93 red SRVs by 
Kiss et al. ( 1999 1 included only the most significant periods 



and found mono-periodic behaviour in 31% of the sample, 
two significant periods for 47%, with a further 13% exhibit- 
ing three dominant periods. By contrast, we have extracted 
all periods above specific S/N thresholds, and include closely 
spaced periods. In order to more easily compare our results 
with earlier studies, we sorted periods for each star into as- 
cending order, and omitted period Pi when Pi+i / Pi < 1.1, 
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since they are likely to represent the same mode (Figure 
right panel). 

Figure [l4| shows that the modal number of periods for 
S/N limits of 10 and 15 is 2, which agrees with the earlier 
studies. Likewise, our distributions are remarkably similar to 
that of Kiss et al.[ ([T999), with the primary difference being 
that, as the S/N constraint is lowered, a small number of 
additional periods have been identified. While the curves 
with S/N ^ 10 are quite similar in both plots, the S/N ^ 5 
and S/N ^ 7.5 curves in the left panel are broader and not 
quite as smooth as those in the right panel, indicating that 
as the S/N threshold is lowered, additional closely spaced 
frequencies (of lower amplitude) are identified. 

If both frequencies were of similar amplitude, simulta- 
neously present, and with long mode lifetimes, we would 
expect to see the effect of beating in the light curves, which 
we do not. Although beyond the scope of the current paper, 
we briefly note that we analysed several stars with close pe- 
riods and well sampled LCs, consisting of both CCD and 
PEP observations. By splitting the datasets (temporally) 
into halves, and using both Fourier and wavelet analysis, 
we conflrmed that the dominant frequencies and amplitudes 
changed between the two segments. The results were con- 
sistent for separate analyses of PEP-only, and CCD-only 
subsets, as well as the combined dataset. Thus it appears 
that amplitude and/or phase modulation is the likely cause 
of the closely spaced periods, consistent with stochastically 
excited oscillations. 

We note that previous analyses of OGLE and MACHO 
photometry of the Magellanic Clouds have used S/N limits 
between 3 and 10 (see, for example, [Soszyriski et al.[[2004| 
Derekas et al.|2006[ ). 



7.5 Period ratios 

Period ratios provide a useful mechanism to compare the 
pulsational behaviour of local SRVs to those in the lower 
metallicity environment of the LMC, since they are inde- 
pendent of distance. Figure [T5| plots a Petersen diagram for 
all multi-periodic stars in our sample (S/N 10). Each red 
point represents the ratio of a longer period to the shorter, 
plotted against the log of the longer period. For compari- 
son, we have plotted the period ratios of LMC MACHO red 
giants ( [Derekas et al.|2006[ grey points), which show a strik- 
ing similarity. The most densely populated region occurs at 
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Figure 15. Petersen diagram plotting log Piong against the ratio 
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hort foi" adjacent periods (red) in galactic M giants 
(S/N ^ 10). Data for the LMC perekas et al.|2006| are plotted 
in grey. 



a period ratio of about 1.3, corresponding to first and sec- 
ond radial overtone modes (P1/P2) in the theoretical model 
of I Wood fc Sebo] ( |1996[ ). The ratio is in excellent agreement 
with the MACHO data, and is one of the most common pe- 



riod ratios found by Soszyhski et al. 1 2004 1 in their analysis 



tifled many more variables, and their Petersen diagram (fig 
11) shows a large number of stars extending to logP ~ 1.1, 
in agreement with our local SRVs. A small period ratio has 
many possible interpretations, including two closely sepa- 
rated oscillation modes, high-overtone radial pulsation, a 
single, slightly varying period, or may be an artefact of am- 
plitude or phase modulation. The former interpretation may 
imply the presence of non-radial oscillation modes ( |Kiss fc| 
Lah 2006 1. However, as noted earlier, our data include multi- 
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Figure 14. The number of periods per star as a percentage of 
the sample, plotted for S/N limits of 5, 7.5, 10, and 15. The left 
panel shows all periods, while in the right panel, closely spaced 
periods have been removed, as explained in the text. 



of period ratios in OGLE observations of the LMC. 

The other significant overlap occurs for period ratios 
close to 1. Although the MACHO data contain few short- 
period (log P < 1.4) stars in this region, the analysis of 
OGLE data by Soszyhski et al. ( 2004 1 appears to have iden- 



ple closely spaced periods that probably represent the same 
oscillation mode, which will result in the small period-ratios 
being over represented in our Petersen diagram, although 
this is somewhat reduced by the use of periods with S/N ^ 
10. 

The area of minimal overlap with LMC data around 
logP ~ 2.4 and period ratios of ~ 1.9 is most likely due to 
our sample selection, which excluded stars with GCVS peri- 
ods longer than 100 d. Although they have no obvious ana- 
logue in the MACHO data, we note that the small group of 
stars with period ratios of ~ 1.9 (and log P ~ 1.6) do appear 



in the larger LMC sample of Soszyhski et al. ( 2004 1. Further- 



more, our period ratios are consistent with those found by 
|Wray, Eyer fc Paczyhski] ( |2004[ ) , who demonstrated that ra- 
tios of ~ 1.0, 1.1, 1.3, and 1.9 are common amongst pulsating 
red giants in the galactic bulge. 



8 CONCLUSION & FUTURE WORK 

A fully autonomous instrument has been designed, built, 
and operated for 5.5 years for a unique, homogenous survey 
of the nearest M giant stars. The main results of this paper 
are the following: 

(i) We have determined the characteristic pulsation peri- 
ods of 247 bright, nearby, M giants, the majority of which 
have not been determined previously, thereby increasing the 
number of such stars by a factor of 4. 

(ii) We provide updated A'-band magnitudes, selected 
from several NIR sources, which we carefully compared to 
determine the most reliable values. We include the effect of 
interstellar extinction and show that circumstellar extinc- 
tion can be ignored. 

(iii) We employ a novel technique to correct residual er- 
rors in our flat-fielding images, and successfully use Optimal 
Image Subtraction for a subset of blended stars, yielding 
scintillation-limited differential photometry to better than 
1% precision. 

(iv) Twenty-three stars were observed concurrently with 
PEP and are in excellent agreement with the CCD data. 
Moreover, 21 stars in our sample had periods determined 
previously, and a comparison with literature values showed 
our periods to be in good agreement. 

(v) All stars in the sample were found to be variable, 
leading us to conclude that the Hipparcos variability flag, 
which was used to select the sample, is a good predictor of 
variability. 

(vi) Many stars in the sample exhibit closely spaced peri- 
ods (period ratios < 1.1), similar to those seen in solar-like 
stars, which probably represent changes in amplitude and/or 
phase of a single oscillation mode. 

(vii) Eighty-seven percent of the stars were found to be 
multi-periodic, with the modal number of periods being 2 
(S/N ^ 10) and 3 (S/N ^ 7.5). 

(viii) The period ratios of local M giants were compared 
to multi-periodic stars in the LMC and Bulge, and were 
found to be in excellent agreement. 

Our motivation for conducting this survey was to signifi- 
cantly increase the sample of nearby M giants to enable un- 
ambiguous identification of the local P-L sequences, and to 
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determine the metallicity dependence of the P-L zero-point. 
These issues will be discussed in Paper II. 
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APPENDIX A: LIGHT CURVES 
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Figure Al. A selection of light curves. CCD and PEP observations are denoted by red(solid) and blue(opcn) points, respectively. 
Hipparcos catalogue numbers are listed in parentheses. 
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Figure A2. Power spectra corresponding to the light curves in Figure [Al] Insets show the spectral window using the same horizontal 
scale as the main plot. A horizontal line shows the mean signal in the frequency range 0.20-0.40 cycle d~^. Peaks corresponding to 
periods listed in Table [3] arc indicated with dotted lines. 
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Figure A3. A selection of light curves. 
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Figure A6. Power spectra of selected stars exhibiting short-period pulsation (top 3 panels). The bottom panel is a light curve of WW 
Pic which exhibits a long secondary period of ~ 370 d and a shorter 36-d period. 



